Spatial quadratic solitons guided by narrow layers of a nonlinear material 



Asia Shapira,^'* Noa Voloch-Bloch/ Boris A. Malomed/ and Ady Arie^ 
^Department of Physical Electronics, School of Electrical Engineering, Tel Aviv University, Tel Aviv 69978, Israel and 

* Corresponding author: asiasapi@post.tau.ac.il 

We report analytical solutions for spatial solitons supported by layers of a quadratically nonlinear 
(X*"^') material embedded into a linear planar waveguide. A full set of symmetric, asymmetric, 
and antisymmetric modes pinned to a symmetric pair of the nonlinear layers is obtained. The 
solutions describe a bifurcation of the subcritical type, which accounts for the transition from the 
symmetric to asymmetric modes. The antisymmetric states (which do not undergo the bifurcation) 
are completely stable (the stability of the solitons pinned to the embedded layers is tested by means 
of numerical simulations). Exact solutions are also found for nonlinear layers embedded into a 
, nonlinear waveguide, including the case when the uniform and localized x^^^ nonlinearities have 

' opposite signs (competing nonlinearities). For the layers embedded into the nonlinear medium, 

stability properties are explained by comparison to the respective cascading limit. 

I. INTRODUCTION 

The use of composite materials and engineered optical media opens ways to new modes of the guided wave propaga- 
■ tion, including self-trapped nonlinear ones, in the form of spatial solitons. Recent reviews summarize results obtained 
O ] along these directions for photonic crystals [l| and quasicrystals 0, quasi-discrete media and nonlinear lattices, 
which feature periodic modulation of the local nonlinearity . 

In ordinary settings, optical solitons are supported by uniform nonlinearities (cubic, quadratic, or saturable), which 
may be combined with a periodic grating, that plays the role of a linear potential of the lattice type 0|, and is 
' necessary for stabilizing solitons against the collapse in the multidimensional geometry On the other hand, 
spatially modulated nonlinearities may themselves induce an effective potential In particular, an interesting issue 
is a possibility to support solitons by localized nonlinearities embedded into linear host media. To introduce the topic, 
we will resort here to a couple of simple models that admit analytical solutions, and thus provide for a direct insight 
^ into specific properties of solitons supported by the localized nonlinearity. 
I— The simplest model of this type was introduced in Ref. 0|, in the form of the nonlinear Schredinger (NLS) equation 
for the wave amplitude u{x, z), with the cubic nonlinearity localized at a single point: 

> iu, + il/2)iP,, + 6ix)\u\''u = 0. (1) 

In terms of the optical transmission, z is the propagation distance and x the transverse coordinate. Obviously. Eq. 
, ([T]) amounts to a linear equation at a; < and a; > 0, supplemented by the jump condition for the derivative at x = 0, 
' which is produced by the integration of the equation around x = 0: 

U:,(a; = +0)-M^(x = -0) = -2|u(a;==0)|^M(2; = 0). (2) 

j ' A family of exact solutions to Eq. ([T]), in the form of peakons, is obvious: 

u{x,z)^{2ky^^exp(ikz'-V2k\x\j, (3) 



O 



- 
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with arbitrary propagation constant fc > 0. This family features degeneracy, as the power (norm) of the solutions 
] does not depend on fc, P = \u(x)\^dx = 1. In particular, the formal application of the Vakhitov-Kolokolov (VK) 
criterion, dP/d^ < 0, which is a necessary stability condition for solitons in self-focusing nonlinear media Q, predicts 
neutral stability of solutions ([3]). In fact, all these degenerate solitons are unstable, collapsing into a singularity or 
decaying, as illustrated by another analytical solution to Eq. ([T|), which explicitly describes the onset of the collapse 
at z ^ -0 |i|: 

V' {x, z) = v^-xo/zcxp i {\x\ - ixaf / {2z) . (4) 

Here > is an arbitrary real constant, and z is negative. The same solution ^ with xq < describes decaying 
solitons at z > 0. The power of this solution is also P = 1, irrespective of the value of xq. 

The solitons may be stabilized if a linear periodic potential is added to Eq. ([T]) The stability is also achieved 
if the single ^-function in Eq. ([iJ is replaced by a symmetric pair, which corresponds to the equation introduced in 
Ref. 

iu, + (1/2)^'.. + [6{x - L/2) + Six + L/2)] = 0. (5) 
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Exact analytical solutions to Eq. ([5]) were found in Ref. [loj for symmetric, antisymmetric and asymmetric localized 
modes. The respective, spontaneous-symmetry-breaking (SSB) bifurcation, which generates asymmetric solutions 
from the symmetric ones, takes place, with the increase of the power, at its critical value Per — (8/9) [1 + (1/3) ln2] w 
0.95. In this model, based on the ideal ^-functions, the bifurcation is degenerate, featuring an ultimately subcritical 
character: branches of the asymmetric solutions go backward as functions of P, up to the state, attained at P = 1, in 
which the entire power is concentrated in an infinitely narrow soliton pinned to one of the two (5-functions. Accordingly, 
these branches are fully unstable. The symmetric modes are stable at P < Pa- and unstable at P > Per, while 
antisymmetric modes are completely unstable, although they do not undergo any bifurcation. 

The degenerate character of the model with the two ideal 5-functions is the price paid for its analytical solvability. 
The degeneracy is lifted if the (5-functions in Eq. ([5]) are approximated by regular expressions: 

6ixTL/2)~^6{xTL/2) = ^cxp{~^^^^^] , (6) 

VTra y a'^ J 

with small regularizing parameter a. The numerical analysis of the regularized model demonstrates that the branches 
of asymmetric states turn forward at some P, which causes the stabilization of the asymmetric solutions past the 
turning points. At a > ag ~ 0.2. the SSB bifurcation becomes supercritical, i.e., the branches of the asymmetric 
solutions go forward immediately after they emerge, being completely stable 

Rather than being represented by a single spot or a symmetric pair, as outlined above, the localized nonlinearity 
may be extended to a periodic lattice of J-functions embedded into the linear medium. The description of stationary 
modes in such a model can be exactly reduced to stationary solutions of the discrete NLS equation (ill . [T2| . which 
has been studied in detail [isj . However, the periodic nonlinearity does not admit asymmetric modes. 

A fundamental role in optics belongs to second-harmonic generating systems based on the quadratic (x^^^) nonlin- 
earity [13, [lB|- In this connection, it is relevant to consider media with one or several narrow x^^^ layers embedded 
into a linear planar waveguide. For the single layer approximated by the respective (5-function, exact solutions in 
the form of peakons, similar to those given by Eq. ([3]), were found in Ref. Unlike solutions ([3]), they are not 

degenerate (the total power depends on the propagation constant), a bigger part of the solution family being stable. 

A new problem, which is considered in the present work, is to find double peakons pinned to a symmetric pair 
of x^^-* delta-functions, cf. Eq. ([5]) for the x^^'' nonlinearity. In this model, we report analytical solutions of all 
the types, viz., symmetric, antisymmetric (as concerns the fundamental- frequency component), and asymmetric ones. 
The corresponding SSB bifurcation is subcritical, but nongenerate (i.e., asymmetric branches eventually turn forward 
as stable ones). Antisymmetric modes do not undergo bifurcations, and turn out to be stable. 

Another new configuration is a nonlinear double layer (alias a dipolar layer) , formed by a fused pair of two narrow 
nonlinear stripes with opposite signs. While it would be very difficult to create such a configuration for the Kerr 
nonlinearity, in x''^' systems it is more feasible, as the sign of the nonlinearity may be changed by reversing the 
orientation of ferroelectric domains accounting for the x^^"* interaction. We consider the double layer described by 
function 5'{x) in front of the x'^' terms. A family of exact solutions for solitons pinned to the double layer is found, but 
they all turn out to be unstable (throughout the paper, the stability is tested via direct simulations of the perturbed 
evolution, in the framework of equations with the ideal (5-functions replaced by their regularized counterparts). 

The most challenging problem is to construct analytical solutions for a soliton pinned to a nonlinear layer embedded 
into a nonlinear waveguide. In this setting, the signs of the localized and uniform x^^^ nonlinearities may be identical 
or opposite. We produce exact solutions of two different types for this case, and test their stability. We also find some 
particular exact symmetric solutions for a pair of nonlinear layers inserted into the nonlinear waveguide. 

The paper is structured as follows. In Section 2, we briefiy recapitulate the peakon solution for the single x''^'' 
layer embedded into the linear medium. In particular, we apply the adiabatic approximation to the description of 
peakons pinned to the layer whose strength slowly varies along the propagation distance. The most essential results 
are reported in Section 3, dealing with the pair of x^^^ layers embedded into the linear medium, including exact 
solutions for asymmetric double peakons. Results for the double layer are presented in Section 4, and the nonlinear 
layer(s) buried into the nonlinear waveguide are considered in Section 5. The paper is concluded by Section 6. 



II. THE x^^^ MONOLAYER EMBEDDED INTO THE LINEAR MEDIUM 



The basic model, with a single narrow channel carrying the x'^'' nonlinearity, can be written as follows: 
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0, (7) 

0, (8) 

where Ai and A2 are local amplitudes of the fundamental-frequency (FF) and second-harmonic (SH) components, ki 
and k2 are the respective wave numbers, k is the nonlinearity coefficients, and Afc the phase mismatch. To reduce 
the number of control parameters, we transform Eqs. ([7]) and ([5|): Ai{X, Z) = {l/2)u{x, z), A2{X, Z) = v{x, z)e^^''^ , 
and rescale the spatial coordinates and coefficients by means of fci, x = X/ki, z = Z/ki, 7 = ^/fci, and Q = 2Afc/fci. 
The resulting normalized equations are 



. dAi 1 
''~dZ ^ 2k[ dX^ 

dA2 1 d^A2 



KS{X/ki)A2Ale-'^''^ 
K5{X/ki)Ale'^^^ 



iuz + {l/2)uxx + l5{x)u*v = 0, 



(9) 
(10) 



The integration of the equations in an infinitesimal vicinity of .t = yields the relations for the jumps of gradients of 
the FF and SH fields at a; = 0, cf. Eq. ©, 



Ux{x = +0) - M^(a- = -0) = -27u*(a; = 0)v(a; = 0), 
Vx{x = +0) - Vx{x ^ -0) ^ ~^u'^{x = 0), 



(11) 
(12) 



while the fields themselves must be continuous across x = 0. 

A family of exact stationary solutions to Eqs. ([9]) and (fTOl) in the form of peakons, similar to solution f[3)). is [li 



±2 [k (4fc + Q)]^/'^ -f-^e'''^ cxp ( -V2k\x\ 



2kj 



exp 



(-V2(4fc + Q)|x| 



Equations Q and f|10p conserve the total power, alias the Manley-Rowe invariant^ 



P 



\u{x)\^dx + A\v{x)\^ 



dx 



For peakon p^. its value is 



Ppcakon = 2^/2 (6fc + Q) Wy/^k + Q 



(13) 



(14) 



(15) 



For Q < 0, this dependence P{k) has a positive slope, dP/dk > 0, at fc > —(5/3, and a negative slope at < fc < — Q/3. 
According to the VK criterion, the peakons should be stable for k > —(3/3, and unstable for < A: < —(5/3. For 
Q > 0, condition dP/dk > holds for all fc > 0, hence the entire peakon family is expected to be stable for (5 > 
|16| . It can be verified by direct simulations that the stability of the peakons precisely complies with the predictions 
of the VK criterion. Note also that expression p5p gives rise to a power threshold: the solitons exist if their total 
power exceeds a minimum value, which vanishes only at (5 = 0, 



Pr. 



P{k = 0) = 2V ^/7^ for g > 0, 
P{k = -Q/3) = 2V^Q/72 for Q < 0. 



(16) 



In the case when the nonlinearity strength slowly varies along the x'^^ layer, i.e.. 7 = 7(2) in Eqs. ^ and (fTO|) . the 
adiabatic approximation may be applied, assuming that the solution given by Eqs. (|13p remains locally valid at each 
value of z, with slowly varying fc(z) which is determined by the conservation of the total power, P = const. Then, 
one immediately finds from Eq. ([TSt 



k[z) = \ 



P^-f^iz) 
72 



6V2 



-\ 2 



-Q 



(17) 



There are limitations on the use of the adiabatic approximation: as seen from Eq. ([T7|) , for (5 < this equation makes 
sense if the expression under the radical is positive, i.e., 7^(2) > 24:\Q\/P'^. For Q > 0, the radical is always real; 
however, in that case there is another constraint, fc > [solution (fT5|) does not make sense for fc < 0]. As follows from 
Eq. (fT7)) . this constraint amounts to 7^(2;) > 8Q/P^. Thus, the adiabatic approximation does not allow the peakon 
to pass points where 7(2;) vanishes. 
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III. THE SYMMETRIC PAIR OF NONLINEAR LAYERS IN THE LINEAR HOST MEDIUM 

A. Formulation 

The modification of Eqs. ^ and ([TO)) for two parallel layers is obvious, cf. Eq. ([5|): 



2 '^xx 



0, 
0. 



(18) 
(19) 



Being a novel model for the x^^'' nonlinearity, it is related to its counterpart itSl) with the cubic nonlinearity through 
the cascad^ng InnU, which corresponds to large positive values of mismatch Q BM- In this limit , one can eliminate 
the SH field, using Eq. ([T^ . v w {j/2Q) S (x — L/2) + 6 {x + L/2) v?, where it is necessary to assume that the ideal 

(5- function is replaced by its regularization 5(x), see Eq. ([B]). Then, the substitution of this approximation into Eq. 
(fT5)) leads to equation 
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1 

2^'^'' ' 2Q 



2 



\uVu = 0. 



(20) 



With the ideal (S-function, Eq. (PO)) does not make sense, as 5'^{x) does not exist. Nevertheless, if 5{x^L/2) are 
taken as smooth approximations, Eq. (|20p is meaningful, being tantamount to the accordingly regularized version of 
Eq. (0). 



B. General analysis 

Stationary localized solutions to Eqs. (|18p and ^9)1 are sought for as 



u{x, z) 



-ikz 



^_gx/2fc(x+L/2)^ at X < -L/2, 
Ai cosh (y2kx^ + A2 sinh (^V2kx^ , at - L/2 < x < +L/2, 
A+e-^(^-^/2), at x>L/2, 



(21) 



^_gV5(3ITQ)(x+L/2)^ at x<-i/2. 



v{x, z) = e-^'''^ < Bi cosh ( y/2 (4fc + Q)x + B2 sinh ( ^2 (4fc + Q)x , at - L/2 < x < +L/2, (22 



B+e-VWQ)(--V2)^ at x>L/2. 

Conditions of the continuity of u and u at a; = ±L/2 make it possible to express the inner amplitudes, Ai,2 and -61,2, 
in terms of the outer ones, A± and B± 



Ai = 
Bi = 

Further, we introduce notations 

Sl 



A++A. 



2 cosh {^^/k/2L 

B++B 



A. 



A.-A^ 



2 sinh ^yk/2L 



, B2 = 



B: - B_ 



2 cosh ^(4fc + Q) /2L 2 sinh ^(4A; + Q) /2i 



(23) 



sinh 



(y2fci) 



S2 



sinh V2(4fc + Q)i 



ci = 1 + coth (\/2fcL) , C2 = 1 + coth (^2 (4fc + g)L) , 



(24) 



71 = 7/V2fc, 72 = 7/ V 2(4/0 + 0) 
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Then, conditions and (|12|) for the jump of the wave functions at points x = ±L/2 give rise to the following 
equations: 

ci^_ - siA+ = 27iA_B_, ciA+ - siA_ = 271A+B+, (25) 
C2B- - S2B+ = 72^2^, C2B+ - S2B_ = 72 . (26) 

Using Eqs. ([2^ . one can eliminate i?± in favor of A±: 

72 (C2A^ + S^Al) 12{C2A^.+S2AI) 

^+ " ^2 _ „2 ' - „2 _ „2 ■ 1^7) 

^2 *2 ^2 ''2 



C. Symmetric, asymmetric, and antisymmetric modes 

Substituting Eqs. ([?7|) into Eqs. (|25p and assuming — 7^ (i.e., that the solution is asymmetric) leads to 
the following expressions for the amplitudes of the asymmetric modes: 

A^A^ = (£i±iH)£i, Al + A' = . (28) 

27172 ^ 2c2 7172 



The asymmetry of the stationary mode is characterized by 



{A+ - A^f _ C1C2 - C1S2 - C2S1 
{A+ + A-f C1C2 - C1S2 + C2S1 ' 



(29) 



On the other hand, Eqs. (|25p and (|26p immediately give rise to the symmetric solutions, 

A+ = A^ = (A±j 



(Cl - Sl) (C2 - S2) 



/symm 



27172 



B+ = (S±)sy„™ = (271)"' (Cl - Sl) . (30) 

Then, Eqs. ([23]) yield, for the symmetric solution, {A2)^y^^^ = {B2)symm = 0' ^^^'^ 

(A \ _ (^±)symm , „ n _ (-^±)symm /„, n 

cosh(^yfc72ij cosh(^V(4fc + Q) /2ij 

The asymmetric solutions emerge, with the increase of L, from the symmetric one as a result of the SSB bifurcation, 
at the point at which the asymmetric solution, as given by Eqs. (|28[) . coincides with its symmetric counterpart pOp . 
Setting, accordingly, e = in Eq. ()29p predicts the location of the bifurcation point: C1C2 = C1S2 + C2S1. Only the 
symmetric solution, given by Eqs. ([30]) and ((3T|) . exists at 

C1C2 < C1S2 + C2S1 , (32) 

where Eq. ([29]) formally yields (A+ — A-)^ < 0, while both the symmetric and asymmetric solutions exist at 

C1C2 > C1S2 + C2S1 , (33) 

when Eq. ^ yields {A+ - A_f > 0. 

As follows from definitions (P^ . condition ([5^ holds in the limit when the two layers merge into one, i — > 0, in 
which case ci « si and C2 ~ S2- On the other hand, at L — ^ cx) definitions (|24p imply ci « C2 ~ 2 and si,2 — > 0, 
i.e., condition ((33| holds in the limit of a very large separation between the layers. Therefore, the SSB indeed occurs, 
with the increase of L, at a particular value of the separation. 

It is also possible to find solutions which are antisymmetric in the FF component, with = — A_ and i?+ = +B^, 
while the SH component remains symmetric. In this case, Eqs. (|23p . and (|26p yield 



^+ — ^ (^±)antis- 



_ / (Ci + Sl) (C2 - S2) 

V 27172 



= i?_ ^ (S±),„,i™ = (271)-' (ci + Sl) . (34) 
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Figure 1: (Color online) Left: the double-barrier structure corresponding, in physical units, to separation L = 200 /^m between 
the two symmetric x^^' layers of width a = 0.92 ^m each. Right: a typical example of stable symmetric solitons. The red 
(taller) and blue (lower) continuous curves depict, respectively, the input for the FF and SH fields, taken as per analytical 
solution given by Eqs. (|2ip - (|27p and (|30p . Chains of dots depict the output produced by simulations of Eqs. p8|) . (|19p . with 
the (5- functions approximated as per Eq. ((6|, over the propagation distance corresponding to z = 100 cm. For this and other 
examples, the scaled wavenumber of the analytical solutions is taken as fc = 1. In the present case, the mismatch is Q = 0. 



V ^Jantisvmm / , • — \ ' V 1/antisyini 

sinh(yfc72ij 

and (^i)a„tisymm = {B2) = 0, cf. Eqs. m and m- 



( ) antisymm 

0172^)' 



cosh (^yj (4fc 



(35) 



D. Numerical results 



The numerical analysis of the dual-layer model aimed to address two issues: the form of the bifurcation diagram 
(subcritical or supercritical), which is implicitly described by the above analytical expressions, and the stability of 
the symmetric, asymmetric, and antisymmetric double peakons. 

We took values of rescaled constants and variables which correspond to fc = 1 and the following typical values of 
physical parameters: at the FF wavelength of 1.0645 /^m, the sample was assumed to be stoichiometric lithium tanta- 
late (SLT) with the e-ee interaction, where two extraordinary waves induce an extraordinary nonlinear polarization. 
Undoing the rescalings which lead to the rescaled notation, it is straightforward to see that rescaled mismatch Q = 1 
corresponds, in physical units, to a very small value, 6 m~^, i.e., our actual results corresponds to the nearly matched 
X*-^^ system. 

Further, the refractive indices and the relevant element of the x*-^-* susceptibility tensor were taken according to Ref. 
[l7| (at 100°C, npF ~ 2.1323, ngn = 2.1999 and ^33 = 12.9 pm/V). In the simulations of the evolution of perturbed 
solutions, the (S-functions were replaced by approximation ([6|) with a ~ a few microns. Finally, taking into regard the 
experimentally measured value of the Kerr coefficient in this material, 712 « 14.6 x 10~^^ cm^/W [ig, it is easy to 
check that, for the physical parameters adopted in this work, the x'"^) nonlinearity is negligible in comparison with 
the x^^^ effects. 

Typical examples of stable symmetric, asymmetric and antisymmetric solitons are displayed in Figs. [Jl [2l and[3l 
respectively. Additional simulations, with strong perturbations added to the input fields (not shown here), demonstrate 
that, as it might be expected, the symmetric solitons are stable before the SSB point and unstable past it. The 
simulations also demonstrate that the antisymmetric solitons are always stable. The stability of asymmetric solitons 
is discussed below. 



7 




Figure 2: (Color online) The same as in Fig. [T] but for a stable asymmetric soliton, with L = 50 ^m, a = 1.06 fj,m, and z — 10 
cm. 




Figure 3: (Color online) The same as in Figs. [T] and [21 but for a stable antisymmetric soliton, with L — 200 /im, a = 1.4 /im, 
and z = 100 cm. 



Making use of Eqs. (p3)) -p9 l) to find the asymmetry e, and calculating the total power of the solutions as per Eq. 
p4p . the bifurcation diagrams were drawn in the plane of (P, e), at different values of mismatch Q. Typical examples 
of the diagrams, displayed in Fig. |4l clearly demonstrate the subcritical character of the SSB bifurcation, similar to 
what was found in model ([S]) with the cubic nonlinearity, which corresponds to the cascading limit of the present 
system. However, unlike that model, the present one, even with the ideal i5-functions, is not degenerate, i.e., the 
branches of the asymmetric solutions go forward after reaching the turning point. 

In accordance with general properties of the symmetry-breaking bifurcations (T9| , one should expect that branches 
of the asymmetric solitons corresponding to de/dP > and de/dP < should be stable and unstable, respectively. 
This expectation was confirmed by direct simulations. In particular, all the asymmetric solitons belonging to the 
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Figure 4: (Color online) The asymmetry parameter of the solitons versus the total power [defined as per Eq. \2Q\ \. for positive, 
zero, and negative mismatch; Q = +1 (a), Q = (b), and Q = — 1 (c). 
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Figure 5: (Color online) The evolution of an unstable asymmetric soliton belonging to the solution branch with de/dP < 
is displayed by means of contour plots of local powers of the FF and SH components in the {z,x) plane. The corresponding 
physical parameters are: a — 1.4 fim, L — 200 fim, and the total propagation distance is z = 100 cm. The rescaled mismatch 
and wavenumber are Q = — 1 and k = 0.26. 



positive-slope branch of the e{P) dependence are stable (Fig. [2] shows an example of such a stable soliton), while a 
typical example of the instability of the branches with the negative slope is displayed in Fig. [5] 
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Figure 6: (Color online) The same as in Figs. [T][3l but for an unstable antisymmetric nearly discontinuous soliton pinned to the 
double nonlinear layer of width a = 7.07 fj,m, shown in the left panel. In this case, the mismatch is Q = 0, the total propagation 
distance is 1 cm, and the initial scaled wavenumber is fc = 1. 



IV. THE MODEL WITH THE DOUBLE NONLINEAR LAYER 



A narrow double layer is formed by two adjacent monolayers with opposite signs of the x^^-* coefficients. Accordingly, 
Eqs. © and ([TU)) are replaced by 

iu, + {l/2)u^^^+"fS'ix)u*v = 0, (36) 
2iv, + il/2)v,., - Qv + {j/2)S'{x)u^ = 0. (37) 



An exact stationary solution to Eqs. (j36p and (|37p can be found in the form of antisymmetric discontinuous solitons, 
cf. expressions (fT3|) for the peakons: 



u{z,x) = ±^27 ^e*''^sgn(x)exp (^-%/2fc|5 
v{z,x) = _^-ie2»fc-sgn(x) exp ^-^2 (4A: + Q)\x^ . (38) 

Note that, unlike the peakons, the amplitudes of these solutions do not depend on fc. 
According to Eq. p4|) . the total power of the discontinuous soliton is 

Pdiscont = + 2 (4A: + Q)"'/'] . (39) 

Expression ((39|) does not give rise to any existence threshold, unlike the peakon solutions [cf. Eq. ([T6|) ]. because 
-Pdiscont(fc — >■ 00) = 0. Obviously, Eq. ((5^ leads to dP/dk < 0, hence the VK criterion predicts that the family of the 
discontinuous solitons is completely unstable. This prediction was confirmed by simulations of the evolution of these 
solitons, see Fig. [6l where strong instability dominates the propagation of the soliton even over a relatively short 
propagation distance, z = I cm. 



V. NONLINEAR LAYERS EMBEDDED INTO A NONLINEAR HOST MEDIUM 

A. A single layer 

If the host medium is itself nonlinear, Eqs. ([9]) and ([T0| are replaced by 

iu, + il/2)u,, + [T + -fSix)]u*v = 0, (40) 
2iv, + {l/2)v,,^Qv + {l/2)[T + jSix)]u' = 0, (41) 
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Figure 7: (Color online) The same as in Fig. O but for the evolution of the unstable double-humped soliton, given by Eqs. 
(|42p -(|44 [l with initial scaled wavenumber k = 1 and Q = —3, pinned to the nonlinear layer of width 7.07 /im. The corresponding 
total propagation distance is z — 100 cm. 



where T and 7 account for the bulk and localized x^^^ nonlinearities, respectively, which may have the same or 
opposite signs, the latter situation corresponding to the competing bulk and localized nonlinearities. Two particular 
exact solutions can be found in this model, following the pattern of the well-known Karamzin-Sukhorukov (KS) 
solutions for x^^"* solitons in the uniform medium j20| . 

The first solution is a straightforward extension of the KS soliton, based on the following ansatz: 

u = e^'^'A sech^ {W {\x\ +^)),v = e^^'^^Bsech^ {W {\x\ + C)) • (42) 

Substituting this ansatz into Eqs. (|40|) . (|4T|) and taking into regard jump conditions ([TT|) . ([T2|) . it is easy to find 
parameters of the exact solution: 

k = -0/3, W = v/-Q/6, A = ±V2B, B = - (g/2r) , (43) 
sinh {2W0 = v/^3Q72 (7/r) , (44) 

With identical signs of 7 and T, Eq. yields ^ > 0, i.e., a single-hump profile of the pinned soliton (|42p . For 

opposite signs of 7 and T, Eq. (|44)) produces ^ < 0, hence the corresponding pinned profile p2)) features a local 
minimum at x = 0, and two maxima at 2; = ±|^|. The soliton is expected to be stable in the former case, and unstable 
in the latter one, when it is pinned by the repelling defect in an unstable position. Direct simulations of the evolution 
of the solitons slightly shifted from the equilibrium positions confirm these expectations. In particular, the instability 
of the double-humped soliton is illustrated by Fig. [T] 

Another type of the pinned soliton can be found in the following form: 

u = e'^'A' [sinh {W {\x\ + , v = e^''''- B' [sinh {W {\x\ + . (45) 

The substitution of ansatz (gS]) into Eqs. (gl]) and ([11]), (HH) produces the following results, cf. Eqs. and 
dUl): 

k = -g/3, W ^ A' = ±V2B', B' = +Q/ {2T) , (46) 

sinh {2WO = - V^3Q/2 (7/r) . (47) 

The solution is nonsingular provided that Eq. (|T7)) yields ^' > 0, which is the case for the opposite signs of F 
and 7. 
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To better understand the meaning of these exact solutions, it is instructive to consider the NLS equation which 
corresponds to the cascading hmit of Eqs. (^0]) and (PT|) : 



+ (l/2)u,,+ f + 75(x) \u\^u = 0, 



(48) 



f = rV (2g) , 7 = 7r/Q + 7V (2g) 5(a;) [recall S{x) is the regularized (5-function ©, cf. Eq. ([201)]. In 

the case of F > 0, i.e.. Q > 0, the exact solution to Eq. which is the counterpart of solution (|^^ . is 



= y^2fc/fe''^^sech (V2fc (|a;| + , sinh (^2%/2fc^^ = 2 (^7/f ) \/2fc, 



(49) 



where fc > is the respective wave number. Solution (|49p is valid for both positive and negative 7, i.e., respectively, 
the attractive and repulsive nonlinear defect in Eq. pS)) . The power of solution is 



P 



+ 00 



dx — I7I ^ 



(2|7|/f)%/2fc-sgn(r)W8(7/f) fc + 1 + 1 



(50) 



It immediately follows from this expression that, for either sign of 7, this soliton family satisfies the VK criterion, 
dP/dk > 0, which suggests that solution (|^^ (with identical signs of 7 and F) is stable in the general case too, when 
the cascading limit does not apply. As mentioned above, this expectation was corroborated by direct simulations (not 
shown here). 

The cascading-limit counterpart of solution corresponds to F < 0, 7 > 0. The respective exact solution to Eq. 
and its power are 



2fc/|f|, 



sinh(V2fc(|x| +0 



, sinh 2V2fc^ =2 7/|r| V2fc, 



(51) 



8 7/f fc + 1- 27/f V2fc + 1 



(52) 



Expression ([5^ does not satisfy the VK criterion, as it yields dP/dk < 0, suggesting an instability of solution ps)) 
in the absence of the cascading limit. Indeed, direct simulations of Eqs. pO)) and (|4T|) confirm that this solution is 
unstable, as shown in Fig. [H 



B. A linear layer embedded into the self-defocusing nonlinear medium (the cascading limit) 



The existence of solution ([51]) to the asymptotic NLS equation (|48|) suggests to consider a similar solution for a 
linear attractive layer (i.e., a usual waveguiding channel) embedded into the medium with the uniform self-defocusing 
cubic nonlinearity (in previous works, solitons pinned by the attractive defect were considered in the NLS equation 
with the self- focusing nonlinearity j2TI.[2^). In terms of the x*^^) system, this solution corresponds to the narrow linear 
channel in the limit of the large negative mismatch. The respective version of the NLS equation is 



iuz + {l/2)uxx + r|upu + ^a6{x)u = 0, 



(53) 



with F < and 70 > 0. In fact, we can fix F = —1 in this case; then, an exact solution to Eq. (|53p for a mode pinned 
by the attractive layer, is [cf. Eqs. ([51]) ] 



2ke'''' 



sinh(y2fc(|x| ,tanh ('\/2fc^) = \/2fc/7 



(54) 



with power P = 70 — v 2fc, which features dP/dk < 0, formally contradicting the VK criterion. However, this criterion 
is not relevant for models with the self-defocusing nonlinearity. Actually, an argument in favor of the stability of 
solution ([54| is the fact that its energy is negative. 



E 



{Ux) 



dx — 7ow^ [x = 0) 



2/3) fc3/2 (70/ V2fc - 1) (270/ V2fc + 1) ' 



(55) 



(recall F = — 1 was fixed), hence the solution has a good chance to represent the ground state of the system. Direct 
simulations of Eq. (|53p corroborate the stability of this pinned mode (not shown here) . 
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Figure 8: (Color online) The same as in Fig. [T] but for the evolution of the unstable soliton, given by Eqs. (|45p - (|47| ) with 
initial wavenumber k — 1 and Q = —3, pinned to the nonlinear layer of width 7.07 fim. The corresponding total propagation 
distance is z = 10 cm. 



C. Two embedded layers 



The extension of the model for the pair of symmetric layers is described by the following version of Eqs. (|40p and 
fflT) [cf. Eqs. (|18p and (fT9| in the case of the linear host medium]: 



iuz + -u^x + i r + 7 



2ivz 



:Vxx - 



r + 7 



5\x- 



L 



L 



0, 



= 0. 



(56) 
(57) 



A particular exact solution to Eqs. (|56p and ([57]) can be found, in the form of a symmetric three-hump structure [with 
a maximum at a; = 0, on the contrary to the double- humped solutioir (|2ip. (|22p . which has a minimum at .t = 0], for 
the competing nonlinearities, F > 0, 7 < 0: 



u{x, z) = Ae 



-ikz 



cosh 



-Q/6i\x\~L) 



at |a-| > L/2, 



cosh ^-Q/6a; , at \x\ < L/2, 



(58) 



z) = Be 



-2ikz 



cosh(^^-Q/6(|a;| - L)^ ^ at \x\>L/2, 
cosh (^/-Q/6x)] ^ at \x\<L/2, 



(59) 



where k, A and B are given by the same expressions (j43|) as in the case of solution (j42|) . 

The highly degenerate nature of this solution is demonstrated by the fact that it satisfies the jump con- 
ditions ([TT|l . at points x = ztL/2 at a single value of the strength of the localized nonlinearity, 7 = 

(recall that both Q and 7 are negative, while F is positive, in the present case). As 



(2F/^/^) sinh i^^^Q/6Lj 

for the stability of the three-hump mode, it may be expected that, with a maximum of the local power set between 
two repulsive nonlinear layers, it is definitely unstable at small distance L between the layers (when they tend to 
merge into a single repulsive element, cf. the instability in Fig. [7]), but it may become stable at larger L, when 
the repulsion from the two separated layers traps the power maximum between them. These expectations have been 
corroborated by direct simulations of Eqs. (|56p and (|57p . In particular, an example of the stable mode, for large 
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Figure 9: (Color online) Left: the profile of the nonlinearity coefficient, for separation L = 200 /iin between two symmetric 
layers of width a — 1.77 ^m each. Right: a stable three-humped solution given by Eqs. (|58p . (|59p . for Q = —3, k = 1, and 
z = 100 cm. 



separation L = 200 /xm, is displayed in Fig. [H] (this solution is stable against large perturbations, which is not shown 
here in detail). 



VI. CONCLUSION 



We have produced several exact solutions for spatial solitons supported by the x'^^' layers embedded into a linear or 
nonlinear planar waveguide. The most fundamental solution describes the full set of families of the symmetric, asym- 
metric, and antisymmetric double-humped modes supported by the symmetric pair of the nonlinear layers inserted 
into the linear medium. The exact solutions describe the subcritical symmetry-breaking bifurcation in this system. 
In addition, particular exact solutions of several types were found for nonlinear stripes running through the nonlinear 
medium, including the case of the competition between the uniform and localized x^^'' nonlinearities. The stability of 
the pinned solitons was tested by means of direct numerical simulations. In the case of the pair of nonlinear stripes 
embedded into the linear waveguide, the character of the (in)stability completely agrees with general principles of 
the bifurcation theory. For the layers embedded into the nonlinear host medium, the results for the stability were 
explained too, with the help of the consideration of the cascading limit. 

The theoretical results reported in this paper call for an experimeirtal realization. As an example, for a small phase 
mismatch Afc = 20 m^^ and a typical nonlinear coefficient of 13 pm/V, the input intensities required to observe the 
soliton in the case of the single embedded layer are about 10^ W/cm^. Such intensities are feasible, as demonstrated 
in Ref. (23| . provided that the necessary nonlinear pattern can be fabricated. 

N.V.B. is an Eshkol Scholar from the Israeli ministry of science, culture and sport. 



[1] J. D. Joannopulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Manipulating light with strongly modulated photonic 
crystals: Molding the Flow of Light (Princeton University Press: Princeton and Oxford, 2008). 

[2] A. N. Poddubny and E. L. Ivchenko, "Photonic quasicrystalline and aperiodic structures", Physica E 42, 1871-1895 (2010). 

[3] F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. Assanto, M. Segev, and Y. Silberberg, "Discrete solitons in optics", 
Phys. Rep. 463, 1-126 (2008). 

[4] Y. V. Kartashov, B. A. Malomed, and L. Torner, "Solitons in nonlinear lattices", Rev. Mod. Phys. 83, 247-306 (2011). 



14 



[5] B. A. Malomed, Z. H. Wang, P. L. Chu, and G. D. Peng, "Multichannel switchable system for spatial solitons", J. Opt. 

Soc. Am. B 16, 1197-1203 (1999). 
[6] B. A. Malomed, D. Mihalache, F. Wise, and L. Torner, "Spatiotemporal optical solitons", J. Optics B; Quant. Semicl. Opt. 

7, R53-R72 (2005). 

[7] B. A. Malomed and M. Ya. Azbel, "Modulational instability of a wave scattered by a nonlinear center", Phys. Rev. B 47, 
10402-10406 (1993). 

[8] M. Vakhitov and A. Kolokolov, "Stationary solutions of the wave equation in a medium with nonlinearity saturation," 

Radiophys. Quantum. Electron. 16, 783-789 (1973). 
[9] N. Dror and B. A. Malomed, "Solitons supported by localized nonlinearities in periodic media", Phys. Rev. A 83, 033828 

(2011). 

[10] T. Mayteevarunyoo, B. A. Malomed, and G. Dong, "Spontaneous symmetry breaking in a nonlinear double-well structure", 

Phys. Rev. A 78, 053601 (2008). 
[11] A. A. Sukhorukov and Y. S. Kivshar, "Spatial optical solitons in nonlinear photonic crystals", Phys. Rev. E 65, 036609 

(2002). 

[12[ Y. Kominis, "Analytical solitary wave solutions of the nonlinear Kronig-Penney model in photonic structures", Phys. Rev. 
E 73, 066619 (2006). 

[13[ P. G. Kevrekidis, The Discrete Nonlinear Schredmger Equation: Mathematical Analysis, Numerical Computations, and 

Physical Perspectives (Springer: Berlin and Heidelberg, 2009). 
[14[ C. Etrich, F. Lederer, B. A. Malomed, T. Peschel, and U. Peschel, "Optical solitons in media with a quadratic nonlinearity", 

Progress in Optics 41, 483-568 (2000). 
[15[ A. V. Buryak, P. Di Trapani, D. V. Skryabin, and S. Trillo, "Optical solitons due to quadratic nonlinearities from basic 

physics to futuristic applications", Phys. Rep. 370, 63-235 (2002). 
[16[ A. A. Sukhorukov, Y. S. Kivshar, and O. Bang, "Two-color nonlinear localized photonic modes", Phys. Rev. E 60, R41-R44 

(1999). 

[17[ I. Dolev, A. Ganany-Padowicz, O. Gayer, A. Arie, J. Mangin and G. Gadret, "Linear and nonlinear optical properties of 

MgO:LiTa03", Applied Physics B 96, 423-432 (2009). 
[18[ S. Ashihara, J. Nishina, T. Shimura, and K. Kuroda, "Femtosecond measurement of nonlinear refraction in periodically 

poled lithium tantalate", in OSA 2002 Nonlinear Guided Wave Conference, paper NLMD41-1. 
[19[ G. looss and D. D. Joseph, 1980, Elementary stability and bifurcation theory (Springer, New York, 1980). 
[20[ Yu. N. Karamzin and A. P. Sukhorukov, "Nonlinear interaction of diffracted light beams in a medium with quadratic 

nonlinearity: mutual focusing of beams and limitation on the efficiency of optical frequency converters", JETP Lett. 11, 

339-342 (1974). 

[21[ X. D. Cao and B. A. Malomed, "Soliton-defect collisions in the nonlinear Schredinger equation", Phys. Lett. A 206, 177-182 
(1995). 

[22[ R. H. Goodman, P. J. Holmes, and M. L Weinstein, "Strong NLS soliton-defect interactions", Physica D 192, 215-248 
(2004). 

[23| S. M. Saltiel, D. N. Neshev, R. Fischer, W. Krolikowski, A. Arie, and Y. S. Kivshar, "Generation of second-harmonic 
conical waves via nonlinear Bragg diffraction", Phys. Rev. Lett. 100, 103902 (2008). 



